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^ ! Abstract 

Q\ . Charm production in deep inelastic scattering has been measured with the ZEUS 

detector at HERA using an integrated luminosity of 120 pb~^. The hadronic 

Q . decay channels D^ — t- Kgn'^, Ajf —^ pKg and A^ —^ Avr^, and their charge 

O . conjugates, were reconstructed. The presence of a neutral strange hadron in 

the final state reduces the combinatorial background and extends the measured 
■ J ■ sensitivity into the low transverse momentum region. The kinematic range is 

^: < pt{D+,A+) < lOGeV, \r]{D+,A+)\ < 1.6, 1.5 < Q^ < lOOO GeV^ and 

0.02 < y < 0.7. Inclusive and differential cross sections for the production of D^ 
mesons are compared to next-to-leading-order QCD predictions. The fraction of 
c quarks hadronising into A^ baryons is extracted. 



The ZEUS Collaboration 

H. Abramowicz^^'"'^, I. Abt'^^, L. Adamczyk^^ , M. Adamus^^, R. Aggarwal^, S. Antonelli^, 
P. Antonioli'^, A. Antonov'^^, M. Arneodo^^, V. Aushev^^'^, Y. Aushev^^'^, O. Bachynska^^, 
A. Bamberger^^, A.N. Barakbaev^^, G. Barbagli^'^, G. Bari'^, F. Barreiro^^, D. Bartsch^, 
M. Basile^, O. Behnke^^ J. Behrl^ U. Behrensl^ L. Bellagamba^, A. Bertolin^s, S. Bhadra^^, 
M. Bindi^, C. Blohmi^, T. Bold^^^ g q Boos^^, M. Borodin^e, K. Borrasl^ D. Boscherini^, 
D. Bot^^, S.K. Boutle^^, I. Brock'^, E. Brownson^^, R. Brugnera^^, N. Briimmer'^^, A. Bruni^, 
G. Bruni^, B. Brzozowska^^, P.J. Bussey^", J.M. Butterworth^\ B. Bylsma^*^, A. Caldwell^^, 
M. Capua^ R. Carlin^s, CD. CatteralP, S. Chekanov\ J. Chwastowski^^-/, J. Ciborowski^^.ai^ 
R. Ciesielskiis.^, L. Cifarelli'^, F. Cindolo^, A. Contin^ A.M. Cooper-Sarkar^^, N. Coppolai^'S 
M. Corradi^, F. Corriveau^", M. Costa"^*, G. D'Agostini^^^ p ^g,! Corso^^, J. de Favereau^^, 
J. del Peso^^, R.K. Dementiev^^ , S. De Pasquale^''', M. Derrick^, R.C.E. Devenish'^^, 

D. Dobur^^, B.A. Dolgoshein^^ , A.T. Doyle^o, V. Drugakov^^, L.S. Durkin^e, S. Dusini^^, 
Y. Eisenberg^^, P.F. Ermolov ^^'^ A. Eskreys^^, S. Fang^^'-', S. Fazio^, J. Ferrando^'^, 
M.I. Ferrero''^ J. FigieP, M. Forrest^o, B. Foster^^, S. Fourletov^^''^'', G. Gach^^, A. GaW^, 

E. Gallo^^ A. Garfagnini^s, A. Geiser^^, I. Gialas^i''', L.K. Gladilin^s, D. Gladkov^^, 

C. Glasman^s, 0. Gogota^s, Yu.A. Golubkov^^, P. Gottlicheri^''^, I. Grabowska-Bo^d^^ 
J. Grebenyuk^^, I. Gregor^^, G. Grigorescu^^ , G. Grzelak^^, C. Gwenlan^'^'"'^, T. Haas^^, 
W. Hain^^, R. Hamatsu"^^, J.C. Hart"^^, H. Hartmann^, G. Hartner^^, E. Hilger^, D. Hochman^'^, 
U. Holm22, R. Hori"^*^, K. Horton^^''^^, A. Hiittmanni^, G. lacobucci^, Z.A. Ibrahimi°, 
Y. Iga^\ R. Ingbir"^*^, M. Ishitsuka*^^ H.-P. Jakob^ F. Januschekl^ M. Jimenezes, T.W. Jones^\ 
M. Jungst^ I. Kadenko^s, B. Kahle^^ B. Kamaluddin i°'t, S. Kananov'^^ T. Kanno^^ 
U. Karshon^^, F. Karstens^^, I.I. Katkov^^'', M. Kaur^, P. Kaur^''^, A. Keramidas^^, 
L.A. Khein33, J.Y. Kim^, D. Kisielewska^^ , S. Kitamura^^''^^ R. Klanner^s, U. Kleini^.'", 
E. KofFeman'^^, D. Kollar^^, P. Kooijman'^^, le. KoroP^, LA. Korzhavina'^^, A. Kotanski^^'^, 
U. Kotzl^ H. Kowalskil^ P. Kulinski^^^ O. Kuprash^^, M. Kuze^^ A. Lee^^, B.B. Levchenko^^, 

A. Levy^^, V. Libovl^ S. Limentani^s, T.Y. Ling^^, M. Lisovyil^ E. Lobodzinskal^ 
W. Lohmann^^, B. Lohr^^, E. Lohrmann^^, J.H. Loizides^^, K.R. Long^'^, A. Longhin^^, 

D. Lontkovskyi^e, Q.Yu. Lukina^^, P. Luzniak^^.aj^ j^ Maeda'^^ S. Magill^ I. Makarenko^^, 
J. Malka^^'''^ R. Mankel^^-'', A. Margotti^, G. Marini'^^, J.F. Martin^°, A. Mastroberardino^, 
T. Matsumoto^^'"', M.C.K. Mattingly^, I.-A. Melzer-Pellmann^^, S. Miglioranzi^^'", F. Mo- 
hamad Idris^*^, V. Monaco'^'^, A. Montanari-*^^ , J.D. Morris^''^, B. Musgrave^, K. Nagano^'^, 
T. Namsoo^^'P, R. Nania^, D. Nicholass^''^, A. Nigro'^^ Y. Ningi^, U. Noor^e, D. Notz^^ 
R.J. Nowak^^ A.E. Nuncio-Quiroz^ B.Y. Oh^°, N. Okazaki^^ K. Oliver^^, K. Olkiewiczi^, 
Yu. Onishchuk^e, Q. Qta^^-"-^, K. Papageorgiu2\ A. Parenti^^ E. Paul^ J.M. Pawlak^^ 

B. Pawlik^^, P. G. Pelfer^^, A. Pellegrino^^ , W. Perlanski^^.aj ^ H. Perrey^^, K. Piotrzkowski^^, 
P. Plucinski^^'"'^, N.S. Pokrovskiy^^ , A. Polini-^, A.S. Proskuryakov^-^ , M. Przybycieii-^^, 
A. RavaF, D.D. Reeder^^ B. Reisert^"^, Z.Ren", J. Repond\ Y.D. RF''*», A. Robertson^^, 
P. RolofF^^ E. Ron29, I. Rubinskyl^ M. Ruspa'^^ R. Sacchi^^ A. Salii^^, U. Samson^ 
G. Sartorelli^, A.A. Savings, D.H. Saxon^o, M. Schioppa^, S. Schlenstedt^^ P. Schleper^^, 
W.B. Schmidke'^^, U. Schneekloth^^ , V. Schonberg^, T. Schorner-Sadenius^^, J. Schwartz^^, 



F. Sciullii\ L.M. Shcheglova^^ , R. Shehzadi^ S. Shimizu'^'''°, I. Singh'^'"', I.O. Skillicorn^o, 

W. Slomiiiskii'^, W.H. Smith5^ V. Sola^^ A. Solano''^ D. Son^^, V. Sosnovtsev^^ , A. Spiridonov^^''?, 
H. Stadie22, L. Stanco^s, A. Stern^^, T.P. Stewart^o, A. Stifutkin^^ P. Stopa^^, S. Suchkov32, 

G. Susinno^ L. Suszyckii^, J. Sztuk22, D. Szuba^s.'-, J. Szubal5•^ A.D. Tapper^s, E. Tassi^'^ 
J. Terron^s, T. Theedt^^ H. Tiecke3^ K. Tokushuku^^'^, O. Tomalak^^, J. Tomaszewska^^'*, 
T. Tsurugai3\ M. Turcato^^, T. Tymieniecka^^'"', C. Uribe-Estrada^^, M. Vazquez^^'", 
A. Verbytskyi^^^ O. Viazlo^^, N.N. Vlasoyi^'", O. Volynets^e, R. Walczak^^, W.A.T. Wan 
Abdullah^o, J.J. Whitmore'^°'''^ J. Whyte^*^, L. Wiggers^^, M. Wing^i, M. Wlasenko^ 
G. Wolfl^ H. Wolfe^^ K. Wronal^ A.G. Yagues-Molmal^ S. Yamada^^ Y. Yamazaki^^-?^ , 
R. Yoshida^, C. Youngman^^, A.F. Zarnecki^^^ l^ Zawiejski^^^ O. Zenaiev^^, W. Zeuner^^'°, 
B.O. Zhautykov25, N. Zhmak^e.^, C. Zhou^o, A. Zichichi^ M. Zolko^^, D.S. Zotkin^s, 
Z. Zulkapli^o 



II 



2 



Argonne National Laboratory, Argonne, Illinois 60439-4815, USA 
Andrews University, Berrien Springs, Michigan 49104-0380, USA 



3 INFN Bologna, Bologna, Italy ^ 

^ University and INFN Bologna, Bologna, Italy ^ 

^ Physikalisches Institut der Universitdt Bonn, Bonn, Germany 

^ H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom 

'^ Panjab University, Department of Physics, Chandigarh, India 

^ Calabria University, Physics Department and INFN, Cosenza, Italy ^ 

^ Institute for Universe and Elementary Particles, Chonnam National University, 

Kwangju, South Korea 
^^ Jabatan Fizik, Universiti Malaya, 50603 Kuala Lumpur, Malaysia 

^^ Nevis Laboratories, Columbia University, Irvington on Hudson, New York 10027, 

USA^ 
^^ The Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sci- 
ences, 

Cracow, Poland 
^^ Faculty of Physics and Applied Computer Science, ACH- University of Science and 

Technology, Cracow, Poland ^ 
^^ Department of Physics, Jagellonian University, Cracow, Poland 

^^ Deutsches Elektronen- Synchrotron DESY, Hamburg, Germany 

^^ Deutsches Elektronen- Synchrotron DESY, Zeuthen, Germany 

^"^ INFN Florence, Florence, Italy 

^^ University and INFN Florence, Florence, Italy 

^^ Fakultdt fiir Physik der Universitdt Freiburg i.Br., Freiburg i.Br., Germany '-" 

'^^ Department of Physics and Astronomy, University of Glasgow, Glasgow, United 

Kingdom 
'^^ Department of Engineering in Management and Finance, Univ. of the Aegean, 

Chios, Greece 
^^ Hamburg University, Institute of Exp. Physics, Hamburg, Germany 

^^ Imperial College London, High Energy Nuclear Physics Group, London, United 

Kingdom ^ 
^^ Institute of Particle and Nuclear Studies, KEK, Tsukuba, Japan ^ 

^^ Institute of Physics and Technology of Ministry of Education and Science of Kaza- 
khstan, Almaty, Kazakhstan 
^^ Institute for Nuclear Research, National Academy of Sciences, and Kiev National 

University, Kiev, Ukraine 

^'^ Kyungpook National University, Center for High Energy Physics, Daegu, South Ko- 
rea "^ 

^^ Institut de Physique Nucleaire, Universite Catholique de Louvain, Louvain-la-Neuve, 

Belgium ^ 
^^ Departamento de Fisica Teorica, Universidad Autonoma de Madrid, Madrid, 

Spain 
^^ Department of Physics, McGill University, Montreal, Quebec, Canada H3A 2T8 
'^^ Meiji Gakuin University, Faculty of General Education, Yokohama, Japan 



III 



^^ nAr\c?r'/~ion hi nrt rttir^ atDfytir^ rt t-^ hto i c"» /^ o I lyi o ■/■/! ■/■/! ( ■fz? /l/*/^o/^/^/j/i l-f n i o ot /i -^ ' 



Moscow Engineering Physics Institute, Moscow, Russia 



33 
34 



Moscow State University, Institute of Nuclear Physics, Moscow, Russia 

Max- Planck- Institut fiir Physik, Miinchen, Germany 

^^ NIKHEF and University of Amsterdam, Amsterdam, Netherlands ^ 

^^ Physics Department, Ohio State University, Columbus, Ohio 43210, USA 

^^ Department of Physics, University of Oxford, Oxford, United Kingdom 

38 INFN Padova, Padova, Italy ^ 

3^ Dipartimento di Fisica dell' Universita and INFN, Padova, Italy ^ 

^^ Department of Physics, Pennsylvania State University, University Park, 

Pennsylvania 16802, USA ^ 
^""^ Polytechnic University, Sagamihara, Japan 

^^ Dipartimento di Fisica, Universita 'La Sapienza' and INFN, Rome, Italy ^ 

^3 Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, United Kingdom ^ 

^^ Raymond and Beverly Sackler Faculty of Exact Sciences, School of Physics, 

Tel Aviv University, Tel Aviv, Israel ^ 
^^ Department of Physics, Tokyo Institute of Technology, Tokyo, Japan ^ 

^^ Department of Physics, University of Tokyo, Tokyo, Japan ^ 

^^ Tokyo Metropolitan University, Department of Physics, Tokyo, Japan ^ 

^8 Universita di Torino and INFN, Torino, Italy 

^^ Universita del Piemonte Orientale, Novara, and INFN, Torino, Italy 

Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 

1A7^ 

Physics and Astronomy Department, University College London, London, United 

Kingdom 

Warsaw University, Institute of Experimental Physics, Warsaw, Poland 

Institute for Nuclear Studies, Warsaw, Poland 

Department of Particle Physics, Weizmann Institute, Rehovot, Israel 
^^ Department of Physics, University of Wisconsin, Madison, Wisconsin 53706, USA 
^^ Department of Physics, York University, Ontario, Canada M3J 1P3 *^ 



50 

51 

52 
53 

54 



IV 



A 
B 
C 

D 

E 
F 



G 



H 



K 



M 



N 



O 



R 



supported by the US Department of Energy 

supported by the ItaHan National Institute for Nuclear Physics (INFN) 
supported by the German Federal Ministry for Education and Research (BMBF), 
under contract Nos. 05 HZ6PDA, 05 HZ6GUA, 05 HZ6VFA and 05 HZ4KHA 
supported by the Science and Technology Facilities Council, UK 
supported by an ERGS grant from the Malaysian government 

supported by the US National Science Foundation. Any opinion, findings and con- 
clusions or recommendations expressed in this material are those of the authors and 

do not necessarily reflect the views of the National Science Foundation, 
supported by the Polish Ministry of Science and Higher Education as a scientific 

project No. DPN/N188/DESY/2009 

supported by the Polish Ministry of Science and Higher Education as a scientific 

project (2009-2010) 

supported by the Japanese Ministry of Education, Culture, Sports, Science and 

Technology (MEXT) and its grants for Scientific Research 

supported by the Korean Ministry of Education and Korea Science and Engineering 

Foundation 

supported by FNRS and its associated funds (HSN and FRIA) and by an Inter- 
University Attraction Poles Programme subsidised by the Belgian Federal Science 

Policy Office 

supported by the Spanish Ministry of Education and Science through funds provided 

by CICYT 

supported by the Natural Sciences and Engineering Research Council of Canada 

(NSERC) 

partially supported by the German Federal Ministry for Education and Research 

(BMBF) 

supported by RE Presidential grant N 41-42.2010.2 for the Leading Scientific Schools 

and by the Russian Ministry of Education and Science through its grant for Scientific 

Research on High Energy Physics 

supported by the Netherlands Foundation for Research on Matter (FOM) 

supported by the Israel Science Foundation 

supported in part by the MINERVA Gesellschaft fiir Forschung GmbH, the Israel 

Science Foundation (grant No. 293/02-11.2) and the US-Israel Binational Science 

Foundation 



V 



" also affiliated with University College London, United Kingdom 
now at University of Salerno, Italy 

^ now at Queen Mary University of London, United Kingdom 

'^ also working at Max Planck Institute, Munich, Germany 

^ also Senior Alexander von Humboldt Research Fellow at Hamburg University, Insti- 
tute of Experimental Physics, Hamburg, Germany 

^ also at Cracow University of Technology, Faculty of Physics, Mathemathics and 

Applied Computer Science, Poland 
3 supported by the research grant No. 1 P03B 04529 (2005-2008) 

^ now at Rockefeller University, New York, NY 10065, USA 

^ now at DESY group FS-CFEL-1 

^ now at Institute of High Energy Physics, Beijing, China 

now at DESY group FEB, Hamburg, Germany 

also at Moscow State University, Russia 
"* now at University of Liverpool, United Kingdom 
" on leave of absence at CERN, Geneva, Switzerland 
° now at CERN, Geneva, Switzerland 
^ now at Goldman Sachs, London, UK 

■? also at Institute of Theoretical and Experimental Physics, Moscow, Russia 
'" also at INP, Cracow, Poland 
' also at FPACS, AGH-UST, Cracow, Poland 
* partially supported by Warsaw University, Poland 
" partially supported by Moscow State University, Russia 
'" also affiliated with DESY, Germany 

"" now at Japan Synchrotron Radiation Research Institute (JASRI), Hyogo, Japan 
^ also at University of Tokyo, Japan 
^ now at Kobe University, Japan 
^ supported by DESY, Germany 
^ deceased 
'"' STFC Advanced Fellow 

nee Korcsak-Gorzo 
°^ This material was based on work supported by the National Science Foundation, 

while working at the Foundation. 
" also at Max Planck Institute, Munich, Germany, Alexander von Humboldt Research 

Award 
"^ now at Nihon Institute of Medical Science, Japan 

"■^ now at SunMelx Co. Ltd., Tokyo, Japan 

"^ now at Osaka University, Osaka, Japan 

"'^ now at University of Bonn, Germany 

"* also at Lodz University, Poland 

"■^ member of Lodz University, Poland 

"•^ now at Lund University, Lund, Sweden 



VI 



al 



also at University of Podlasie, Siedlce, Poland 



VII 



1 Introduction 

Heavy-quark production in ep interactions in deep inelastic scattering (DIS) is dominated 
by the Boson Gluon Fusion (BGF) process. Heavy-quark production provides a two-fold 
test of perturbative quantum chromodynamics (pQCD); a study of the BGF process and 
the higher-order corrections to it, and an independent check of the validity of the gluon 
density in the proton extracted from the inclusive DIS data. Of the two heavy quarks 
whose production is accessible at HERA, c and b, the latter is strongly suppressed due to 
its smaller electric charge and larger mass. 

The production of charm via the identification of D and D* mesons in DIS has been exten- 
sively studied at HERA in the kinematic range 1 < Q^ < 1000 GeV^ Pt{D, D*) > 1.5 GeV 
[IHH], where Q^ is the negative squared four-momentum exchange at the electron vertex 
and pt is the transverse momentum. The results are consistent with the calculations of 
pQCD. The fragmentation fraction /(c — )■ A^) has been measured by the ZEUS collabo- 
ration in the photoproduction regime [9]. The obtained fragmentation fraction is larger 
than but consistent within uncertainties with the average from e~^e~ collisions |10] . 

In this paper, a charm quark in the final state was identified by the presence of a charmed 
hadron. The production of D^ mesons and A+ baryons was studied using the decay^j 
D~^ — )■ KgTi^, A^ — !■ pKg and A+ — ?■ Avr^. These decay channels were chosen since the 
presence of a neutral strange hadron in the final state significantly reduces the combina- 
torial background. Measurements of D~^ and A+ cross sections provide information about 
both c-quark production and its fragmentation. 

With respect to previous studies, in this analysis the kinematic region of the measurement 
is extended to very low transverse momenta of the produced charmed hadrons. No explicit 
cut on the transverse momenta of the reconstructed charmed hadrons was applied. This 
is particularly relevant at low Q^, where charm quarks are predominantly produced with 
low transverse momentum. In addition, A+ production was studied for the first time at 
HERA in DIS. From a comparison of the D~^ and A^ cross sections, the fragmentation 
fraction /(c — )■ A+) is extracted. 



2 Experimental set-up 

The analysis was performed with data taken from 1996 to 2000 corresponding to a lumi- 
nosity of 120.4 ± 2.4 pb~^. The sample consists of 38.6 pb~^ of e~^p data collected at a 



^ The charge conjugated modes arc iniphed throughout this paper. 



centre-of-mass energy of 300 GeV and of 65.1 pb collected at 318 GeV, plus 16.7 pb of 
e~p data collected at 318 GeVo 

A detailed description of the ZEUS detector can be found elsewhere [TT]. A brief outline 
of the components most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [12], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consisted 
of 72 cylindrical drift chamber layers, organised in 9 superlayers covering the polar-anglqj 
region 15° < 9 < 164°. The transverse-momentum resolution for full-length tracks was 
(^{Pt)/pt = 0.0058pr © 0.0065 © 0.0014/pt, with pt in GeV. 

To estimate the ionisation energy loss per unit length, dE/dx, of particles in the CTD [13] , 
the truncated mean of the anode-wire pulse heights was calculated, which removes the 
10% lowest and at least the 30% highest pulses depending on the number of saturated 
hits. The measured dE/dx values were corrected by normalising to the measured average 
dE/dx for tracks around the region of minimum ionisation for pions with momentum p 
satisfying 0.3 < p < 0.4 GeV. Henceforth, dE/dx is quoted in units of minimum 
ionising particles (mips). 

The high-resolution uranium-scintillator calorimeter (CAL) [H] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
was subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAG). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a{E)/E = 0.18/vE for electrons 
and a{E)/E = 0.35//E for hadrons, with E in GeV. 

The position of the scattered electron at the CAL was determined by combining in- 
formation from the CAL and, where available, the small-angle rear tracking detector 
(SRTD) [15] and the hadron-electron separator (HES) [T6] . 

The luminosity was measured from the rate of the bremsstrahlung process ep — )■ e'jp, 
where the photon was measured in a lead-scintillator calorimeter [17] placed in the HERA 
tunnel a.t Z = —107 m. 



^ Hereafter, both electrons and positrons are referred to as electrons, unless explicitly stated otherwise. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity is 
defined as 77 = —In (tan 2), where the polar angle, 9, is measured with respect to the proton beam 
direction. The azimuthal angle in the X-Y plane is called 0. 



3 Theoretical predictions 

The next-to-leading-order (NLO) QCD predictions for the cc production cross sections 
were obtained using the HVQDIS program [18] based on the fixed-flavour-number scheme 
(FFNS). In this scheme, only light quarks (n, d and s) and gluons are included in the 
proton parton density functions (PDFs) which obey the DGLAP equations [19], and the 
cc pair is produced via the BGF mechanism [20] with NLO corrections [21]. The presence 
of different large scales, Q, px and the mass of the c quark, rric, can spoil the convergence 
of the perturbative series because the neglected terms of orders higher than al (where as 
is the strong coupling constant) contain log(Q^/m^) factors which can become large. The 
FFNS variant of the ZEUS-S NLO QCD fit [22] to structure function data was used as the 
parametrisation of the proton PDFs. In this fit, as{Mz) was set to 0.118 and the mass of 
the charm quark was set to 1.5 GeV; the same mass was used in the HVQDIS calculation. 
The renormalisation and factorisation scales were set to jJiR = jji-f = \/Q'^ + 4m^. The 
charm fragmentation to the D^ meson was modelled using the Peterson function [23] with 
the Peterson parameter, e, set to 0.079 [21]. For the hadronisation fraction, /(c — )■ -D^), 
the value 0.2161q;o29 ^as used [7]. 

The HVQDIS predictions for the production of D^ mesons are affected by the theoret- 
ical uncertainties listed below. The uncertainty on the total cross section is given in 
parentheses: 

• the ZEUS PDF uncertainties were propagated from the experimental uncertainties of 
the fitted data (+5:25); 



the charm quark mass was changed consistently in the PDF fit and in HVQDIS by 

-15.2%\ 
-13.5%^' 



±0.15GeV(+i^-2^"^ 



the renormalisation scale was varied by a factor 2 (]']i2g%)' 

the factorisation scale was changed by a factor 2 independently of the renormalisation 

scale (_2i,7%J) 

the e parameter of the Peterson fragmentation function was changed to 0.01 and 

0.1 [2H[25]. This modification affects the shapes of the pt, Q^ and x distributions 



4 Monte Carlo models 

The detector acceptance was modelled using the Rapgap 3.00 [2B] Monte Carlo (MC) 
program, interfaced with Heracles 4.6.1 [27] in order to incorporate first-order elec- 
troweak corrections. The generated events were passed through a full simulation of the 



detector, using Geant 3.13 [28], and finally processed and selected in the same way as 
the data. 

The MC was used to simulate events containing charm quarks produced in the BGF pro- 
cess. The Rapgap generator used leading-order matrix elements with leading-logarithmic 
parton showers. The CTEQ5L [25] PDFs were used for the proton. The charm-quark 
mass was set to 1.5 GeV. Charm fragmentation was simulated using the Lund string 
model [30]. The D^ and A+ hadrons originating from beauty decays were accounted for 
by including a Rapgap 6-quark sample where the 6-quark mass was set to 4.75 GeV. An 
additional sample where charm was produced by the process eg — t- eg was generated and 
was used to study the model dependence of the simulation. For this process, the charm 
quark was treated as a part of the structure of the photon. The processes gg — )■ ec and 
qq — )■ cc were not included because their contribution estimated using the Rapgap MC 
was found to be less than 1% in the studied kinematic range. 

In general, the MC gives a reasonable description of the data for DIS and D~'"-meson 
variables when compared at detector level. To improve the description further, Rapgap 
was reweighted to reproduce the pt{D^) distribution observed in the data. The same 
weights used for D^ mesons were also applied to D'^ and A+ hadrons. 



5 Kinematic reconstruction and event selection 

A three-level trigger system was used to select events online [TT|[3T]. At the third level, 
an electron with an energy greater than 4 GeV and a position outside a box of 24 x 12 cm^ 
centred around the beampipe on the face of the rear calorimeter was required by a fully 
inclusive DIS trigger which had a high acceptance for Q^ ^ I GeV^. However, this trigger 
was heavily prescaled and the equivalent luminosity is 17pb~^. 

Additionally, events above Q^ ~ 20 GeV^ were selected by a medium-Q^ trigger. The only 
difference to the inclusive DIS trigger is that the position of the scattered electron on the 
RCAL face had to lie outside a circle centred around the beampipe of radius between 25 
and 35 cm, depending on the running period. 

The fraction of the electron energy transferred to the proton in its rest frame, y, as well 
as the kinematic variables Q^ and Bjorken x, were reconstructed offline using the electron 
method [32] (denoted by the subscript e), which uses the energy and angle of the scattered 
electron. The inelasticity y was also obtained using the Jacquet-Blondel (JB) method [53] . 
The double angle (DA) method [32], which relies on the angles of the scattered electron 
and the hadronic-energy flow, was used as a systematic check. 

The following requirements were imposed offline: 



38 < 6 < 65 GeV, where 6 = ^ Ei{l — cos 9i) and Ei and 9i are the energy and the 
polar angle of the i'^'^ energy-flow object (EFO) [31] reconstructed from charged tracks, 
as measured in the CTD, and energy clusters measured in the CAL. The sum i runs 
over allEFOs [35]; 

E^ > 10 GeV, where E^ is the energy of the scattered electron identified using a 
neural-network algorithm 



• E'conc < 5 GeV, where Econe is the calorimeter energy measured in a cone around the 
electron position that was not assigned to the electron cluster. The cone was defined 

by Rcone < 0.8 with i?eone = ^(^0)^ + (Ar/)2; 

• a match between the tracking and the calorimeter information for electrons well within 
the CTD acceptance, 17° <6e< 149°. For 6e outside this region, the cut 5 > 44 GeV 
was imposed; 

• for events with the scattered electron reconstructed within the SRTD acceptance, the 
impact position of the electron on the face of the RCAL had to be outside the region 
26 X 14 cm^ centred on X = Y = 0. If the electron position was reconstructed without 
using SRTD information, a box cut of 26 x 20 cm^ was imposed; 

• 1.5<Ql< 1000 GeV^; 

• yjB > 0.02 and y^ < 0.7; 

• a primary vertex position in the range | ^vertex I < 50 cm. 

This analysis used charged tracks measured in the CTD that were assigned either to the 
primary or to a secondary vertex. The tracks were required to have transverse momenta 
Pt > 0.15 GeV and pseudorapidity in the laboratory frame |?7| < 1.75, restricting the 
study to a region where the CTD track acceptance and resolution were high. Candidates 
for long-lived neutral strange hadrons decaying to two charged particles were identified by 
selecting pairs of oppositely charged tracks, fitted to a displaced secondary vertex. The 
events were required to have at least one such candidate. 



6 Strange-particle reconstruction 

The Kg mesons were identified by their charged decay mode. Kg — )■ tt+tt^. Both tracks 
were assigned the mass of the charged pion and the invariant mass, M{n^7T~), of each 
track pair was calculated. Additional requirements to select Kg were imposed: 

• M{e^e~) > 50MeV, where the electron mass was assigned to each track, to eliminate 
tracks from photon conversions; 



• M{p7t) > 1121 MeV, where the proton mass was assigned to the track with higher 
momentum, to ehminate A contamination in the Kg signal; 

• cos 6'xY > 0.98, where 9xy is defined as the angle between the momentum vector of 
the Kg candidate and the vector defined by the primary interaction vertex and the 
Kg decay vertex in the X-Y plane; 

• 483 < M(7r+7r-) < 513 MeV; 
. \r,iKl)\ < 1.6. 

The A candidates were reconstructed by their charged decay mode to pvr". The track with 
the larger momentum was assigned the mass of the proton, while the other was assigned 
the mass of the charged pion, as the decay proton always has a larger momentum than 
the pion, provided the A momentum is greater than 0.3 GeV. Additional requirements to 
select A were imposed: 

• M(e+e-) > 50 MeV; 

• M(7r+7r^) < 483 MeV, where the charged pion mass was assigned to both tracks, to 
remove Kg contamination in the A signal; 

• cos^xy > 0.98; 

• 1112 < M{pn) < 1121 MeV; 

• \ri{A)\ < 1.6. 

Figure [T] shows the invariant-mass spectra of Kg, A and A candidates. Distributions of 
the reconstructed proper lifetime for these particles based on the same data sample as 
analysed in this paper were found to be satisfactory [38] . 



7 Reconstruction of charmed hadrons 

The production oi D^ and A+ hadrons was measured in the range of transverse momentum 
< pt(-D"^,A+) < 10 GeV and pseudorapidity |?7(Z}+, A+)| < 1.6. Strange-hadron candi- 
dates were combined with a further track measured in the CTD which was assigned to 
the primary interaction vertex. The combinatorial background was significantly reduced 
by requiring pT{D~^)/E!p'^^ > 0.1 and pt(A^)/-E'^'^^° > 0.12, where the transverse en- 
ergy iJy-^^^ was evaluated as El^^^ = Xlie >io°(-^« ^^^^«)- '^^^ ^^^ Tuns over all energy 
deposits in the CAL with a polar angle 6 above 10°. The details of the reconstruction of 
the three different decay channels are given in the next subsections. 



7.1 Reconstruction of the decay D^ -^ I^^tt"'" 

The D^ mesons were reconstructed from the decay channel D^ — )■ A'°7r+. In each event, 
D^ candidates were formed from combinations of K'g candidates reconstructed as de- 
scribed in Section with further tracks assumed to be pions. The pion candidates were 
required to have pt{t^^)/E!^^^^ > 0.04. Only pion candidates with dE/dx < 1.5 mips 
were considered. Further reduction of the combinatorial background was achieved by 
cutting on the angle between the pion in the D~^ rest frame and the D^ flight direction, 
6*{'ir~^). Different cuts depending on pt{D~^) were used to ensure optimal background 
suppression: 

• cosr(7r+) < 0.9 for 0.0 < Pt{D+) < 1.5 GeV; 

• cos^*(7r+) <0.8 for 1.5 < prP+) < 3.0 GeV; 

• cosr(7r+) < 0.6 for 3.0 < PtP+) < 10.0 GeV. 

The Kg7r~^ invariant-mass distribution was fitted with the sum of contributions from the 
signal, the non-resonant background and a reflection caused by Df — )■ KgK~^ decays. 
The signal was described by a Gaussian function defined as: 

g{a, Mo; m) = ^=^ exp \ , '\ (1) 

Vivrcr no- 



where Mq and a are the resonance mass and width, respectively. For the background a 
sum of Chebyshev polynomials up to the second order was used: 

b{A, B, C- y{m)) = A-{l + B-y + C- {2y' - 1)), (2) 

where y{m) = (2m - mmax - "^min) / ("^max - m^m) and mmax("^min) = 2.1(1.6) GeV is 
the upper (lower) limit of the fitted range. 

The mass distribution of the reflection r^m) caused by the decay D'^ -^ KgK^ — )■ 
7i^7i~K~^ was obtained from D'^ combinations in the Monte Carlo at detector level 
matched to the same decay at generator level. The normalisation of the reflection with 
respect to the Gaussian signal assumed for D^ — )■ Kgii^ decays is based on previously 
measured fragmentation fractions / [7] and branching ratios B [32] (see also Table [T]) and 
the detector acceptances for both decay channels. For this purpose, the invariant mass 
distribution of the reflection was normalised to unity and then multiplied by the expected 
ratio of Dt to D~^ mesons: 



/(c -^ Dt) ■ B{Dt -^ K^sK^ -^ 7T+n-K+) AjDj 
f{c -^ D+) ■ B{D+ -^ Kl'K+ -^ 7r+7r-7r+) A{D^ 



R = ^;; .^^i;^;::. X. -.11. ■ ^r)^ = 0-44 ± 0.10, (3) 



where A{D'^) and Ai^D'^) are the reconstruction acceptances for D^ and D~^ mesons, 
respectively, as obtained from the Monte Carlo. The resulting fitting function is given by: 

F{A, B, C, D, a, Mq; m) = b{A, B, C; y{m)) + D ■ [r{m) + g{a, Mq; m)], (4) 

where the parameters A, B, C, D, a and Mq were determined by the fit. 

Figure [2] shows the invariant mass spectrum for the D~^ candidates after the reflection 
was subtracted using the fit, resulting in a 20% reduction in the number of D~^ mesons. 
A clear signal is visible. The fit yielded a D~^ mass of 1872 ± 4 MeV, in agreement with 
the PDG value [39]. The width of the signal was 19.0 ±3.1 MeV, refiecting the detector 
resolution. The number of D~^ mesons yielded by the fit was N{D~^) = 691 ± 107. 

In order to extract the D"'"-meson yields in bins oi p'^{D~^), ri{D~^), Q^ and x, the signals 
in all analysis bins of a given quantity were fitted simultaneously, fixing the ratios of the 
widths in the bins to the Monte Carlo prediction. All other parameters including the 
masses were left free for all bins in the simultaneous fit. 

The signal in the region < pt{D^) < 1.5 GeV that was not accessible in previous 
measurements is shown in Fig. [31 



7.2 Reconstruction of the decay A+ — )- pKg 

The A+ baryons were reconstructed from the decay channel A^ — )■ pKg. In each event, A^ 
candidates were formed from combinations of Kg candidates reconstructed as described 
in Section [6] with proton candidates. The proton-candidate selection used the energy-loss 
measurement in the CTD. Tracks fitted to the primary vertex with more than 40 hits 
were considered. The proton band was parametrised separately for positive and negative 
tracks from an examination of dE/dx as a function of the momentum |3D]. The proton 
selection was checked by studying proton-candidate tracks from A decays. To remove the 
region where the proton band completely overlaps the pion band, the proton momentum 
was required to be less than 1.5 GeV and a cut on dE/dx > 1.2 mips was applied. Due 
to the proton selection described above, refiections from D~^ — )> Kg-K^ and Df -^ KgK~^ 
decays are suppressed. 

As a result of the cut on the proton momentum, there is no acceptance for A^ baryons 
at very high pr(A^). Hence the measurement of the cross section for this decay channel 
was restricted to the region < pr(A^) < 6 GeV. 

Figure [H shows the M{pKg) distribution for the A+ candidates. A clear signal is seen at 
the nominal value of the A^ mass [39]. The mass distribution was fitted to the sum of 
a Gaussian function describing the signal and the function defined in Eq. ([2]) to describe 



the non-resonant background. The number of reconstructed A^ baryons yielded by the 
fit was iV(A+) = 79 ± 25. 

7.3 Reconstruction of the decay A+ —^ Att"'" 

The A+ baryons were also reconstructed from the decay channel A^ — >• Avr^. In each event, 
A+ candidates were formed from combinations of A candidates as described in Section El 
with further tracks assumed to be pions. The pion candidates were required to have 
Pxij^^) / E^^^ > 0.05. Only pion candidates with dE/dx < 1.5 mips were considered. To 
suppress combinatorial background further, the cut cos6'*(7r"'") < 0.8 was imposed, where 
6*{7r^) is the angle between the pion in the A+ rest frame and the A+ fiight direction. 

Figure |5] shows the M(A7r) distribution for the A+ candidates. Wrong-charge combinations 
in the data sample, normalised to the right-charge combinations in the region outside the 
peak, are also shown. For wrong-charge combinations, the sum of the charges of the 
proton from the A candidate and the further track is equal to zero. The data were fitted 
to the sum of a Gaussian function describing the signal and the background function 
defined in Eq. ([2]). The number of reconstructed A+ baryons obtained from the fit was 
A^(A+) = 84 ± 34. 

The signal-to-background ratio for both studied A+ decay channels is similar. Figure [6] 
shows the invariant-mass spectrum containing both A^ -^ pKg and A+ — ;■ Avr"*" candi- 
dates. The fit yielded A^(A+) = 146 ± 33 candidates. This combined peak was not used 
to extract any cross sections or fragmentation fractions. 



8 Cross sections and acceptance corrections 

For a given observable, Y, the differential cross section in a bin i was determined using 

daj _ iV,(D+) 
dY Ai-C-B-AYi' 

where Ni{D^) is the number of reconstructed D+ mesons in bin i having size AFj. The re- 
construction acceptance, Ai, takes into account migrations, efficiencies and QED radiative 
effects for the i^^ bin, C is the integrated luminosity and B is the branching ratio [3H] for 
the decay channel used in the reconstruction (see Table [1]). The total visible production 
cross sections were determined using 

A-C-B ' 



where N{D~^, A^) and A were determined for the whole kinematic range of the measure- 
ment. All acceptances were obtained from the Monte Carlo. 

The 6-quark contribution, predicted by the MC simulation, was subtracted from all mea- 
sured cross sections. The Rapgap prediction for beauty production was multiplied by 
two, in agreement with a previous ZEUS measurement of beauty production in DIS |41j . 
The subtraction of the 6-quark contribution reduced the measured cross sections by 2 — 3% 
for the D~^ and about 1% for the A^. 

There is no sizeable acceptance for charmed hadrons in the transverse-momentum range 
< pt{D~^,A.'^) < 0.5 GeV. Hence an extrapolation using the reference Monte Carlo 
was performed when the cross sections were extracted. For example, the extrapolation 
accounts for 6% of the D'^ production in the full kinematic range of the measurement and 
for 11% of the D^ production in the restricted range < pt{D~^) < 1.5 GeV. 



9 Systematic uncertainties 

The systematic uncertainties of the measured cross sections and fragmentation fractions 
were determined by changing the analysis procedure and repeating all calculations. In the 
measurement of the differential and total cross sections, the following groups of systematic 
uncertainty sources were considered. The effects on the total cross sections are shown in 
parentheses {D+; A+ -^ pK'^; A+ -^ A7r+): 

• {6i} event and DIS selection {'ts%'-i -2%' -4%)- '^^^ following cut variations were applied 
to data and MC simultaneously: 

— the cut on i/jb was changed to yjB > 0.03; 

— the cut on the scattered electron energy E^ was changed to E'g > 11 GeV; 

— the cuts on 6 were changed by +2 GeV; 

— the cut on | ^vertex | was changed to | Z^crtcx \ < 45 cm; 

— additionally, a box cut of 26 x 14 cm^ was used for all electron candidates without 
an SRTD requirement; 

• {^2} Q^ and X reconstruction (<1%; —3%; —6%). The DA method was used for the 
reconstruction of Q^ and x instead of the electron method; 

• {^3} energy scale (±2%; ^^y^; ^4%). To account for the uncertainty of the absolute 
CAL energy scale, the energy of the scattered electron was raised and lowered by 1% 
and E^^^^° was raised and lowered by 2%. These variations were only applied to the 

MC; 

• {64} model dependence of the acceptance corrections: 
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— the process eg — t- eg was included in the Rapgap MC sample (+5%; +3%; +9%); 

— the MC samples were not reweighted in priD^^ , D+ , A.+) (-17%; -6%; -21%); 

• {5{\ uncertainty of the beauty subtraction (Ig^; ±1%; <1%). This was determined 
by varying the subtracted 6-quark contributions by a factor 2; 

• {5q} uncertainty of the signal extraction procedure (^9% ; ^5% ; ^8% )• 

— the fit was repeated changing the invariant mass window of 1.6 — 2.1 GeV by 
±50 MeV on both sides for D^ — )■ K^tt^ decays. Similarly, the considered invariant 
mass region of 2.0 — 2.5 GeV was changed by ±50 MeV for A+ — )■ pK^ decays and 
by ±30 MeV for the channel A+ -^ A7r+; 

— the choice of the background function was assigned an uncertainty of ±5%. This 
value was estimated by comparing the fit results obtained using different choices for 
the background function, such as polynominals of different orders or exponential 
functions; 

— for differential cross sections, the assumed Gaussian width ratios were varied by 
±10%; 

• {5j} uncertainty in the luminosity measurement of ±2.0%. 

The following uncertainty was considered only for the decays D^ — )■ K^-k^ and A+ — )■ K^p: 

• {^s} K^ reconstruction (±2%; ±1%; — ). Since the MC signal had a narrower width 
than observed in the data, the invariant-mass window for the K^ candidate selection 
was reduced to 0.486 < M(7r+7r-) < 0.510 GeV in the MC only. 



The following source of uncertainty was considered only for the decay D^ — ?■ K^ti^: 

• {^9} uncertainty of the reflection subtraction(±5%; — ; — ). The normalisation of the 
D^ reflection was changed by the uncertainty of R (see Eq. ([3])) due to the uncertainties 
of the fragmentation fractions and branching ratios used in the calculation. 

The following source of uncertainty was considered only for the decay A+ — )■ K^p: 

• {<^io} proton reconstruction (— ; —14%; — ). The following checks were performed: 

— the number of hits required for the proton candidates was lowered to 32; 

— the uncertainty of the dE/dx simulation for low-momentum protons was evaluated 
changing the parametrisation of the proton band |40] : 

— the cut on the energy loss was lowered to dE/dx > 1.15 mips. 

The following source of uncertainty was considered only for the decay A+ — )■ Avr^: 

• {^ii} A reconstruction (— ; — ; ±4%). Since the MC signals had a narrower width than 
observed in the data, the invariant-mass window for the A candidate selection was 
reduced to 1.113 < M{p7r) < 1.120 GeV in the MC only. 
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Contributions from the different systematic uncertainties were calculated and added in 
quadrature separately for positive and negative variations. These estimates were made in 
each bin in which the differential cross sections were measured. Uncertainties due to those 
on the luminosity measurement and branching ratios were only included in the measured 
D^ and A^ total cross sections. For differential cross sections, these uncertainties are not 
included. 

As an additional check, the dE/dx efficiency for pions and protons was verified directly 
in the data using K^ and A decays. For the D^ — )■ K^-k^ decay channel, the effect of 
the dE/dx cut on the pion candidate tracks was very small and the result changed only 
marginally when the cut was released. 

The average cross sections obtained from the two different running periods {^/s = 300 
and 318 GeV) are expressed in terms of cross sections at ^/s = 318 GeV. This involves a 
typical correction of +1% determined using HVQDIS. 



10 Results 

Charm hadron cross sections were measured using the reconstructed D^ and A+ signals 
(see Section E]) in the kinematic range < pt{D+,A^) < 10 GeV, |?7(L'+,A+)| < 1.6, 
1.5 < g2 < 1000 GeV^ and 0.02 <y< 0.7. 

In addition to the statistical and systematic uncertainties, a third set of uncertainties is 
quoted for the measured cross sections and charm fragmentation fractions, due to the 
propagation of the relevant branching-ratio uncertainties (Tabled]). 

10.1 D'^ cross sections 

The following total visible cross section for D^ mesons was measured: 

(7(L>+) = 25.7 ±4.1 (stat.) +|^ (syst.) ± 0.8 (br.) nb. 

The corresponding prediction from HVQDIS is a{D^) = 12.7 ^^4^ nb. The measured and 
predicted cross sections are in agreement to better than two standard deviations. 

To allow a direct comparison to a recent measurement of D~^ production by the ZEUS 
collaboration using a lifetime tag [8], the cross section was extracted for the kinematic 
region defined by 1.5 < PTiD+) < 15 GeV, |r/(D+)| < 1.6, 5.0 < Q"^ < 1000 GeV^ and 
0.02 < y < 0.7. The measurements using different decay channels and different techniques 
were found to be consistent. 
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The differential cross sections as functions of p\{D^), ri{D~^), x and Q^ are shown in 
Fig. [7] and given in Table |2J The cross sections in Q^ and x fall by about three orders of 
magnitude, while the cross section in p'^{D^) falls by about two orders of magnitude in 
the measured region. There is no significant dependence of the cross section on ri{D^). 
The HVQDIS predictions describe the shape of all measured differential cross sections 
reasonably well. The differential cross section in p^i^D^) is compared to a previous ZEUS 
result [7j for p\{D'^) > 9 GeV^. The two measurements are in good agreement. 

10.2 A+ cross sections and fragmentation fractions 

The following A^ cross sections were measured: 

• using the decay channel A^ — )■ pKg in the restricted range < pr(A^) < 6 GeV: 

a(A+) = 14.9 ±4.9 (stat.) tH (syst.) ±3.9 (br.) nb; 

• using the decay channel A+ — t- Att^: 

a(A+) = 14.0 ± 5.8 (stat.) tH (syst.) ± 3.7 (br.) nb. 

To compare and combine both measurements, the value obtained for the decay channel 
A^ — 7- pKg was multiplied by 1.01 ± 0.01 to extrapolate to the full kinematic region 
considered in this paper. The cross sections obtained using different decay channels are 
in good agreement. To extract the A^ fragmentation fraction, the measurements were 
combined taking into account all systematic uncertainties and their correlations: 

acombincd(A+) = 14.7 ±3.8 (stat.) tli (syst.) ±3.9 (br.) nb. 

The uncertainty of the branching ratio was treated as partially correlated since both 
branching ratios, B{A'^ — )■ pKg) and B{A^ — )■ Att"*"), were measured relative to the decay 
mode A+ ^ pK-n^ [35] • 

The fragmentation fraction f{c^ A+) can be calculated using the D'^ cross section: 

/(c^A+) = ^^-/(c^D+). (5) 

(j[D+) 

In a previous ZEUS publication [7] /(c — ?■ D'^) was defined as: 

where a^{D^), cr'^(Z}°) and a^{D'^) are the cross sections for Pt{D) > 3 GeV. The factor 
1.14 takes into account the production of charm-strange baryons [7]. For D^ and D^ 
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mesons the equivalent cross sections (as described elsewhere [9]) were used. Combining 
Eqs. (E]) and dSJ yields: 

Since the cross sections cr(-D^) and cr(A^) were measured down to pt{D~^ , A'^) = GeV, 
no treatment of the different transverse momentum distributions for D~^ and A^ hadrons 
was necessary. The measured value: 

/(c -^ A+) = 0.117 ± 0.033 (stat.) to^oH (syst.) ± 0.027 (br.), 

is compared to previous measurements in Table [31 The result is consistent with a previous 
ZEUS measurement in the photoproduction regime [5] and with the e~^e~ average value. 



11 Conclusions 

Open-charm production in ep collisions at HERA has been measured in deep inelastic 
scattering using three decay channels. The presence of a neutral strange hadron in the 
final state allowed the measurement to be extended to very low transverse momenta of 
the reconstructed charmed hadrons. The total visible and differential cross sections for 
D~^ production are in reasonable agreement with NLO QCD predictions. The measured 
D~^ cross sections are consistent with previous ZEUS results. The fragmentation fraction 
f{c^ A^) has been measured for the first time at HERA in deep inelastic scattering. The 
result obtained from a combination of two decay channels is consistent with a previous 
measurement performed in the photoproduction regime and with the average e~^e~ value. 
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Decay mode 


Branching ratio [%] 


D+ -^ K°s7r+ -^ n+n-n+ 
D+ -^ K+Kl -^ K+n+n- 

A+ -)■ A7r+ -)■ p-K^n^ 


1.00 ±0.03 
1.03 ±0.06 
0.80 ±0.21 

0.68 ±0.18 



Table 1: Branching ratios of the charmed hadron decays j[3^ 



pI{D+) bin 


da/dp%{D+) 


^stat 


^syst 


( GeV') 


(nb/ GeV') 


(nb/ GeV') 


(nb/ GeV') 


0, 2.25 


7.1 


±2.1 


±1.3 -1.1 


2.25, 4.41 


3.3 


±0.9 


±0.4 -0.3 


4.41, 9.0 


0.80 


±0.22 


±0.17 -0.16 


9.0, 100.0 


0.026 


±0.007 


±0.004 -0.006 


r/(D+) bin 


da/dr]{D+) 


^stat 


^syst 




(nb) 


(nb) 


(nb) 


-1.6, -0.5 


7.5 


±1.9 


±1.1 -1.5 


-0.5, 0.5 


6.8 


±1.6 


±0.9 -1.8 


0.5, 1.6 


10.3 


±2.6 


±1.9 -1.9 


g^bin 


da/dQ"^ 


^stat 


^syst 


( GeV^) 


(nb/ GeV') 


(nb/ GeV') 


(nb/ GeV') 


1.5, 5.0 


4.0 


±1.3 


±1.0 -0.5 


5.0, 40.0 


0.33 


±0.06 


±0.03 -0.06 


40.0, 1000.0 


0.0013 


±0.0004 


±0.0003 -0.0002 


X bin 


do jdx 


^stat 


^syst 




(nb) 


(nb) 


(nb) 


0.000021, 0.0004 


43000 


±12000 


±9000 -8000 


0.0004, 0.0016 


7300 


±1400 


±800 -1400 


0.0016, 0.1 


19.2 


±5.7 


±2.8 -3.7 



Table 2: Measured D~^ cross sections as a function ofp\{D^), ri{D^), Q^ andx for 
1.5 < g2 < 1000 GeV"^, 0.02 <y< 0.7, < Pt{,D+) < 10 GeV and |r/(D+)| < 1.6. 
The statistical and systematic uncertainties are shown separately. The cross sec- 
tions have further uncertainties of 3% from the D^ — )■ K^ii:^ — )■ ti^t\:~tx^ branching 
ratio, and 2% from the uncertainty in the luminosity measurement. 





f{c ^ A+) 


ZEUS (DIS) 


0.117 ±0.033 (stat.) t^^^ (syst.) ± 0.027 (br.) 


ZEUS (7p) [g 


0.144 ±0.022 (stat.) t','^ (syst.) t',^, (br.) 


combined e^e~ data 


0.076 ± 0.007 (stat. © syst.) +°:°27 (br.) 



Table 3: The fraction of c quarks hadronising to a A+ baryon, f{c — )■ A^) 
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Figure 1: Mass distributions of the secondary vertex candidates in the (a) Kg, (h) 
A and (c) A samples. The statistical uncertainties are in general smaller than the 
point size. For illustration the data have been fitted using the sum of a '^modified" 
Gaussian function ^9^ and a linear background. 
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Figure 2: The M{Kgn~^) distribution (dots) for D^ candidates. The reflection 
caused by the decay D^ — )■ K^K^ has been subtracted as described in the text. 
The solid curve represents a fit to the sum of a Gaussian signal and a background 
function, while the background contribution alone is given by the dashed curve. The 
dotted histogram shows the reflection scaled as described in the text with an offset 
of 680 to position it at the bottom of the figure. 
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Figure 3: The M{K^tt'^) distribution (dots) for D^ candidates in the region 
< pt{D^) < 1.5 GeV . The reflection caused by the decay Df — )■ KgK~^ has been 
subtracted as described in the text. The solid curve represents a fit to the sum of 
a Gaussian signal and a background function, while the background contribution 
alone is given by the dashed curve. 
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Figure 4: The M{pKg) distribution (dots) for A+ candidates in the region < 
Pt(A^) < 6 GeV . The solid curve represents a fit to the sum of a Gaussian signal 
and a background function, while the background contribution alone is given by the 
dashed curve. 
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Figure 5: The M(A7r"^) distribution (dots) for A^ candidates. The solid curve 
represents a fit to the sum of a Gaussian signal and a background function, while the 
background contribution alone is given by the dashed curve. The dotted histogram 
shows the distribution of wrong-charge combinations (see text). 
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Figure 6: The invariant mass distribution (dots) for A^ — )■ pK^ and A+ -^ Ayr"'" 
candidates. The solid curve represents a fit to the sum of a Gaussian signal and a 
background function, while the background contribution alone is given by the dashed 
curve. 
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Figure 7: Differential D^ cross sections as a function of (a) p'^{D^), (b) r]{D^), 
(c) Q^ and (d) x compared to the NLO QCD calculation of HVQDIS. The measured 
cross sections are shown as dots and the triangle represents a previous ZEUS result. 
The X-axis in (a) is broken. The inner error bars show the statistical uncertainties 
and the outer error bars show the statistical and systematic uncertainties added in 
quadrature. The band shows the estimated theoretical uncertainty of the HVQDIS 
calculation. 
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